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Abstract: We study the distance-redshift relation in a universe filled with 'walls' of 
pressure-less dust separated by under dense regions. We show that as long as the density 
contrast of the walls is small, or the diameter of the under dense regions is much smaller 
than the Hubble scale, the distance-redshift relation remains close to what is obtained 
in a Friedmann universe. However, when arbitrary density contrasts are allowed, every 
prescribed distance-redshift relation can be reproduced with such models. 
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1. Introduction 

Since more than a decade, cosmology research is facing the dark energy problem: the 
present Universe seems to be in an accelerating phase. This conclusion was first drawn from 
measurements of the distance-redshift relation from type la Supernovae (SNIa) |l], |2| and 
is confirmed by many other datasets, from the cosmic microwave background Q to baryon 
acoustic oscillations and other aspects of large scale structure. Until very recently the 
measurements inferring the existence of dark energy rely mainly on the distance-redshift 
relation which is valid in a Friedmann Universe New independent measurements of, 
e.g. the expansion rate H{z) are now being performed see e.g. Q. Hence this situation 
is changing, so that we shall soon know both, dj^(z) and H{z) with good accuracy. The 
general opinion is that fluctuations on large scales are small so that they can be treated 
with linear perturbation theory and linear perturbations average out in the mean over many 
directions and large scales, and therefore fluctuations are not relevant for the determination 
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of quantities like cIa(z) and H(z). This expectation has been confirmed by perturbative 
calculations. Within linear perturbation theory, the fluctuations of the distance-redshift 
relation for redshift z > 0.2 is on the level of a few percent Q. 

However, perturbations on smaller scales can become very large, density fluctuations 
e.g. in galaxies are Sp/p ~ P g a,i/p io 8 . Since the relation between metric perturbations, 
or more precisely the Christoffel symbols, and density fluctuations is non-linear, it is not 
evident that large amplitude, non-linear, small scale density fluctuations cannot add up to 
affect the distance-redshift relation on large scales. 

To study the real problem one would need to analyse light rays passing through a 
realistic Universe with high density fluctuations. So far, this has been done only within 
Newtonian N-body simulations, see e.g. However, it is well known that Newtonian 
gravity misses the terms which are relevant for the back reaction problem ||, hence a full, 
non-linear relativistic treatment is needed. Since this is very difficult, so far mainly toy 
models which mimic reality to a certain extent have been studied. 

The present work inscribes in this framework. Instead of considering spherically sym- 
metric solutions of general relativity (GR), the so called Lemaitre ||-Tolman JlCfl-Bondi 
(LTB) models, for recent reviews see we study a Universe containing high density 
walls. We shall consider infinitely extended parallel walls. The considered model is a sub- 



case of the Szekeres solution [12[. Light propagation in general Szekeres model has been 
studied recently ]l3|, 14]. This is of course a gross over-simplification, but we know that 
galaxies tend to be aligned in filaments and photons coming to us from a far away super- 
nova, might experience a geometry similar to the one of such a symmetric wall universe. 
The weakest point of our toy model is that all the walls are parallel while we expect a 
typical photon to traverse filaments which are aligned in different directions. We shall take 
this into account to some extent by studying photons coming in from different directions 
with respect to the walls. 

Such walls have been studied in the past but only perturbatively. Since we know 
that the effects are small within linear perturbation theory, we cannot trust higher order 
perturbation theory if it predicts large deviations from the Friedmann distance-redshift 
relation. For this reasons we analyse exact, fully relativistic wall-universes in this work. 

In the next section we present the wall metric and the Einstein equations. We also 
study the conditions on the parameters which have to be satisfied so that no singularity 
apart from the Big Bang is present in the backward light cone of the observer. In section || 
we present the results for the distance-redshift relation for 'realistic' walls and for a wall 
universe which mimics the observed relation. In section ^] we conclude. 

2. Wall Universes 

In this section we study universes containing only pressure-less matter (dust) and which are 
symmetric under translations and rotations in a plane which we call the y-plane. They have 
the same number of symmetries as LTB models and can be solved analytically, see ]l6| . 
The metric is of the form 

ds 2 = -dt 2 + a 2 (t,x)dx 2 + b 2 (t,x)(dyl + dy\) . (2.1) 
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Note that the only difference to the LTB geometry is that our symmetrical 2d manifolds 
are planes, dy\ + dy\ = dr 2 + r 2 dcj) 2 while those of LTB are spheres, dJ7 2 = d6 2 + sin 2 9d<p 2 . 
We denote the spatial coordinates by x = (x,y 1,1/2) in order to reserve the letter z for 
the redshift. In the following a prime denotes a derivative w.r.t. x while a dot denotes 
derivative w.r.t. t. The Einstein equations for this geometry and for pure dust matter 
yield [0, m 



d t [-)=d t E = 0, (2.2) 



a 



; ., i'b'\ 2 M(x) , . 

b ~ - = 2 — r 1 > ( 2 - 3 ) 



ay 6 

M' = Ai:Gpb 2 b' = ^Gpb 2 aE(x) . (2.4) 
In Eq. ( |2.2|) we have introduced the time-independent function 

E(x) = b'/a (2.5) 



and Eq. ( |2.3| ) defines M(x) which is also time-independent. In LTB models M/G can be 
interpreted as mass density (Note that in the LTB case a term 6/(2(7) has to be added 
to M which is a consequence of the curvature of the 2-sphere. For more details see fli~7f| .) ; 
and (M' /G)r 2 dr is the mass in a shell of thickness dr. However as the mass in an infinite 
plane is not well defined, this interpretation is not meaningful in the planar case. In our 
case it is therefore not unreasonable that M may become negative even though a, b and p 
are supposed to be positive at all times. 

From the matter conservation equation we also obtain dt(pb 2 a) = 0, which, on the 
other hand, is a consequence of Eq. (|2.4| ). 

2.1 The solutions 

Eq. (|2.3| ) can we rewritten as 



j , 2 = 2MW +£(l)2) {%6) 



with parametric solutions [12, [0] 



M 2M n 

for E + : b = (coshr? - 1) = -p- sinh 2 (r?/2) , (2.7) 

M 

t = — (sinh 77 — 77) +t B (x), forM>0; (2.8) 
E 6 

b = -|1 (cosh 77 + 1) = (sinh 2 (r7/2) + l) , (2.9) 
M 

t = -—. r(sinh 77 + 77) + t B (x), for M < ; (2.10) 
E 6 

b = \E\{t - t B {x)) forM = 0; (2.11) 
3 _ \ 2 /3 



for£ = 0: b = \^-V2M(t - t B {x))\ , for M > , (2.12) 
b = b = const. , for M = . (2.13) 



3 



Note that for E = Eq. ( |2,q ) implies that M > 0. This equation also implies 



b > 



2M 



at all times, in all cases. 

The function is arbitrary; it is called the 'bang time'. For M > 0, at t = tg, i.e 

rj = 0, we have 6 = which represents the Big Bang singularity. Positions with M < 
have no Big Bang singularity but a 'bounce' at t = t B - We shall simplify below to the 
case t B = 0, i.e., uniform bang time. Note that we have chosen expanding solutions. From 
these we can obtain the collapsing solutions simply by changing the sign of t. Since in the 
Einstein equations only b appears they are invariant under t — >• —t. 

Of course the {t =const.} hypersurfaces are not parallel to the {rj =const.} hyper- 
surfaces, but their position depends on x. For fixed position x, Eqs. @P and 
correspond to Friedmann solutions with curvature K = —E 2 < and M = AirGpb 3 /3. 
Note that unlike in the Friedmann case, wall solutions with M < need not be unphysical. 

The parametric representation with r/ is chosen in order to express the solutions in 
terms of elementary functions, but it is of course not necessary. For example, for M > 0, 
setting 

/ t \ 2 / 3 

r(t,x) = -E^—) and 
S(t) = (-3T) _1 sinli 2 Q [sinh-id]" 1 (6(-r) 3/2 )^ 



we obtain 



b(t,x) 



Note that in the definition of S, [sinh —id] 1 denotes the inverse of the function in brackets, 
and id is the identity function, id(x) = x. One can check that S solves the differential 



equation [18| 



i(S + rS') 2 + 3r- i = 0, 



(2.14) 



1 /3 

with initial condition S(0) = (|) . Note that this is the only regular solution, i.e solution 
with 5"(0) oo. This expression will be useful in Section |3.3[ 



The function a(x,rj) can be obtained from Eq. (2.5). For example for M > we find 



E- 1 



db 

Ox 



for E / : 
for E = : 



a = — f — - | sinh 2 (— 
E\E 2 \2 



coth 



t' B + 



M 
£3 



(t ~ t B ) 2/3 

M l/3 6 l/3 



E 



+ 



5 x 6V3 



(sinh rj— rj) , (2.15) 



(2.16) 



(The suffix t in db/dx indicates that we have to interpret b as functions of (t, x), not (a;, rj), 
in this derivative.) Even if E = 0, Eq. (2.4) implies that < M'/E < oo, so that the 
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r.h.s. of Eq. ( [2.16 ) is well defined. Below, we shall choose the x-coordinate such that 



M'/E =constant. 

Note that M[x) and E{x) can pass through zero so that in general different solutions 
from above have to be glued together at the boundary of their validity. We have checked 
that this gluing process can be performed in a smooth way and does not induce singularities 
in the scale factor b. However, for M — >• the scale factor a — > oo. Nevertheless, we believe 
this to be a coordinate singularity, since, as we have checked, both, the Kretschmann scalar, 
K = R a j3^ v R a ^^ u an d the scalar curvature remain finite for M — > 0. In our examples below 
we shall have M > throughout and therefore we do not encounter this problem. However, 
when computing a from Eq. (|2.5| ), one has to be careful to use the result ( 2.15| ) and take 



the limit E — > for fixed t, hence also 77 — ^ 0. One cannot use ( |2.12| ) and ( [2.5D , since 
for E = we have M' = so that Eq. ( |2.5| ) is identically satisfied and cannot be used to 
obtain a(t, x). 

2.2 Singularities 

Singularities can occur when a, b or p become either infinite or zero. To have no singularities 
(apart from the Big Bang) which occurs at t = tB, hence b = 0, in the past light cone of 
every possible observer we might be interested in, we must demand that all singularities 
lie in the future. In more precise models, when one specifies the observer location, one can 
relax this condition to the one that no singularity lies within the background lightcone of 
the specific observer. 

In general, the question of singularities depends on the choice of the functions M(x) 
and E(x). From our solutions it is clear that b behaves monotonically as a function of time 
for fixed x. This is to be expected since no clustering goes on in the directions y\ and 1/2 
described by this scale factor. Since we are interested in an expanding b, a singularity is 



present when the the scale factor a of the x-direction tends to zero. From Eq. ( 2.15 ) we 
infer that for ts = 0, a = implies 

2 (M/E 2 )' cosh(r ? /2) 



E (M/E 3 )' sinh 3 (?7/2) 



(sinhr/ — 77) > 4/3 . 



It is easy to verify that the right hand side is an even positive function with minimum 4/3 
at rj = 0. Hence there is a singularity at some finite value of r] if the l.h.s. ever becomes 
larger than 4/3 or, equivalently, if 

E ' M/E 3 _ cosh (?? /2) 

E(M/E 3 )'- 2sinh 3 (r / /2) (Smhr? V) ^ 1/6 

for some value of x. 

We now consider a simple ansatz motivated by the perturbative analysis presented in 



Ref. [16]. We choose 

2 



9t, 

and 



M(x) = ^ (1 + eh{x)) (2.17) 



AirGpb 2 a = — = -t^ = const. (2.18) 
E 3 
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so that 

h! 

E = ej. (2.19) 

In full generality M'/E = fix) could be an arbitrary positive function of x. But we 
can always make a coordinate transformation to x{x) determined by 

dx 1 

dx 6irGpb 2 atQ ' 

so that with respect to the new coordinate M'/E =constant. Hence we just fix the co- 
ordinate x (up to a constant shift) by this choice. In addition, we have chosen uniform 
bang time, ts(x) = 0. This is a true restriction. With this we have reduced the three 
free functions of x to one, h(x) which defines the density profile. Furthermore, we have 
introduced the parameter e such that for e = we reproduce the matter dominated Fried- 
mann solution. We may also require |/i(x)| < 1 so that e indicates the amplitude of the 
perturbations. We do this in one of the examples below. 

The above requirement for a singularity at some time t ^ now reduces to MM" < 
M /2 /3. (Strictly our derivation applies only for M' ^ 0. For M' oc E = 0, one sees 



directly from Eq. Q2.16 ) that M" oc E' < is the necessary and sufficient condition for 
a = at some time t > ts-) We have found that most interesting mass profiles satisfy 
this condition for some values of x and therefore have singularities at some time in some 
places. This is not surprising but actually expected from gravitational collapse. However, 
when over densities become very high and we approach the collapse, pressure forces and 
heating become important and our simple pressure-less dust model for matter no longer 
holds. In order to be able to stay within the present framework, we therefore demand 
that such singularities be in the future and not in the past for the density profiles under 
consideration. 

Let us consider ELS db first example 

h{x) = cos(kx) . 

Then the condition for the existence of a singularity (at t ^ 0) becomes 

— (ecos(/cx) + e 2 cos 2 (A;x)) < (e 2 /3) sin 2 (A;x) , 

which is always satisfied for some values of x, irrespective of k and e. A similar behavior is 
expected whenever h is not a convex function, but a function representing several under- 
and over-densities cannot be convex. 

However, this is not so important for our considerations. As we have said, the require- 
ment of singularities to be absent is mainly a technical one and it is actually sufficient not 
to have a singularity in the past. 



Using the above expression for a (for M > 0) and the ansatz ( 2.17 , 2.19 ) for M and E, 
we find that a = is equivalent to 

(1 + eh)h" 1 1 1 cosh(77/2) . , , . . . 

' ' (sinhr/ - rj) - 1 > -1/3 . (2.20) 



eh*-3(l + eh)h» 3l_ 5(T ^_ 2 s inh 3 (^/2) 
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Interestingly, in extremal positions of h, with b! = 0, the l.h.s. of the above expression 
is exactly —1/3. This comes from the fact that for this case r] = V t and we have to 
replace the condition that there is no singularity before some given time to by a{t) > for 
t < to using expression (2.16) for a(t). If h!" = when h! = (as in our example) one can 
show that in the positions where h has a maximum, hence b! = and h" < 0, 1 + eh > 0, 
singularities occur first. Furthermore, when 1 + eh > and h" < 0, the denominator of 
the l.h.s. of Eq. (2.2C) is larger than 1 and hence the l.h.s. becomes larger than —1/3. 
Therefore, there exists a finite value r] s (x) where Eq. ( 2.20| ) is satisfied and a(x,rj s (x)) = 0. 
If, on the contrary, 1 + eh > and h" > the l.h.s. of Eq. ( [2.20 ) is smaller than —1/3. For 
positions in the vicinity of an extremum this implies that if the extremum is a minimum 
of h, the position x does not encounter a singularity in the future while positions close to 
maxima do. 

Let us study in more detail the request that the second singularity (not the big bang 
one) lies in the future, t > to- Using the expression ( |2.8| ) for t, we can rewrite the condition 
a(x, 7] s ) = as 

(1 + eh)h" _ 9 cosh(?7 s /2) tlt(x,r] s )e 3 h' 3 



eh' 2 - 3(1 + eh)h" 4 sinh 3 (7? s /2) (1 + eh) 
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The condition t(x,rj s ) > to, for h' < which we shall consider hence r/ s < for t(x,rj s ) > 0, 
then becomes 



(1 + eh) 
e 3 h' 3 



1 + 



(1 + eh)h" 



ih' 2 



9tl 



< 



cosh(7?s/2) 
sinh 3 (ry s /2) 



3(l + eh)h" 

This equation for r] s (x) can only be solved numerically. However, often we realize that the 
l.h.s. is smallest at small \h'\ i.e. for small values of |2£(:r)|. Hence singularities will develop 
first in positions with small \h'\. This requires also small \rj a \ so that we may develop the 
scale factor a and t in r] s . The above inequality then leads to power law relations and 
inserting the above expression for E = (3/2)M'tg yields the constraint 



1 + 



(3t2) 7 /3 3 l/3 2 2/3 



80 



M' 2 \ 

MVS 



QM"M 1 / 3 

1 - 7^(^) 2/3 (6ecos(fcr) (1 + e cos(A;x)) 1/3 



> 0. 



+ 



(sin (kx)Y 



> 0. (2.21) 



(l + ecos(A:a;)) 2/3 ' 

The first inequality is general while for the second inequality we have chosen h = cos(kx). 
In Fig. [l] we plot the constraint for this case together with the condition to use the limiting 
solution for E = 0, ( |2.16| ), (which is not necessary for our analysis) in the e-A plane, where 
A denotes the wavelength of the perturbation A = 2ir/k. 



3. The distance redshift relation in a wall universe 

3.1 Generalities 
3.1.1 Redshift 

We now consider a photon emitted from a source at some position and time (t s , x s ) arriving 
in our telescope at position and time (to,xo). We denote the matter 4-velocity field, hence 
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<g 1000 




Figure 1: The region above the red line has singularities in the future only. While the blue 
line describes the condition to use the limiting solution for E = 0, ( 2.16| ). This can be used when 
tE 3 /M <§; 1, where with "<C" we mean at least two orders of magnitude smaller. The two black 
lines describe the physical parameters e = 9.5 x 10~ 6 and A = 80 Mpc. The green line is the 
Hubble scale Hq 1 . With physical parameters we mean an amplitude as determined by WMAP |H| 
observations and a wavelength agrees with the size of the largest observed voids |2(J which is about 
40-90 Mpc. More precisely we find e requiring that at early time there is only a single density 
fluctuation in each Hubble distance. This leads, at first order, to S — 87r 2 e/15, and the matter 
density fluctuation at early times, 6 = 5 x 10 -5 can be inferred from WMAP observations. For 
more details see |l5f. 



the 4- velocity of source and observer by u(t, x) and the photon 4- velocity by n. The redshift 
of the source, z is then given by 

l + z= y ) ' !\ a . 3.1 

We consider a co-moving source and observer, hence u = dt and normalize the affine 
parameter of the photon, s, such that n°(so) = 1. The redshift then reduces to 

l + z = n°\ s (3.2) 

for our geometry with goo = — 1 and goi = 0. From the geodesic equation for the photon 
we infer that its momenta in y\- and ^-direction are simply redshifted so that 

J\ = \?n} = b 2< ^- = const. and J<i = b 2 n 2 = b 2 -^- = const. (3.3) 
ds ds 



hence 



n \ 1 
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From the geodesic equation for n we can now derive the evolution of the redshift: 
dz_ dn» _ 2 y Jj + Jl fb b>\ 

Here we have used a = b'/E to eliminate the scale factor a. Note also that the prime and 
the dot in the above equation denote partial derivatives while d/ds is a total derivative 
along the path of the photon. 

3.1.2 Distance 

The evolution of the distance to the source is given by the Sachs focussing equation pl| , 

d 2 D 



ds 2 



\a\ 2 +TZ)D. (3.6) 



D is the angular diameter distance to the source, a is the complex scalar shear of the light 
bundle which we define below and 

n = ^R^n u = AnGT^rfn" = 4vrG(l + z) 2 (p + P) . (3.7) 

Here P is the pressure in the direction of the photon. The important point is that this 
quantity is non-negative for any energy momentum tensor which satisfies the dominant 
energy condition p > P in all directions, hence also for a cosmological constant where we 
have TZ = 0. In terms of the affine parameter of the photon, the growth of the angular 
diameter distance to the source is not accelerated. If the dominant energy condition is 
satisfied D(s) is always a concave function. Furthermore, clustering which leads to the 
production of non-vanishing shear is only increasing the deceleration of D as function of 
the affine parameter s. But of course we do not measure this function but D{z) which can 
behave very differently. 

The complex shear of the light ray bundle is defined as follows (appendix 7.7.3. in [p^]): 
We consider two spatial orthonormal vectors e\ and e 2 which are normal to both, u and n 
at the observer and are parallel transported along n, such that V n e a = for a = 1, 2. The 
vectors ei,e2 are a basis of the so called 'screen'. Note that we do not require that u be 
parallel transported along n, hence ei,e2 are in general not normal to u elsewhere than at 
the observer, where we have given their initial conditions. The complex shear is defined by 

a = ^g(e, V e n) , e = ei + ie 2 . (3.8) 

In order to compute the shear we must know n not only along the photon geodesic itself 
but we must determine its derivatives in directions normal to n. We shall directly use 
the transport equations p^ j. For a vorticity free ray bundle (which is the case here) with 

b 

9 + 9 2 + a\ + a 2 = -TZ, (3.9) 

a x + 20CTi = -Re (J 7 ) , (3.10) 

& 2 + 29a 2 = Im {?) , (3.11) 



expansion rate 9 = \n^ >1 these are 
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where a\ = Re (a), 02 = Im(<r), and T = \R a ^ v ( OL (^ri il n v . To determine the shear a 
we need to know the initial conditions for the differential equations ( |3.9| ) to (3.11). It 



is possible to determine the behavior of the shear and the expansion of the light near the 
vertex [23|. Choosing the affine parameter of the photon to vanish at the observer position, 



so = 0, these are 



a{s) = - S -T + O(s 2 ), (3.12) 



1(i-1k sA +0{s 3 ). (3.13) 



J-"o and 1Zq are the values of T and 1Z at the observer position. The light bundle expansion 
9 diverges at the observer position, but we can consider an initial condition not exactly at 
the observer. This choice can affect the numerical precision. After determining 1Z, e and T 
for a given geometry and photon direction, we can solve the system ( |3.9| ) to ( |3.11 ) together 
with the Sachs focusing equation (|3.6| ) numerically. 

3.2 'Realistic' walls 

We want to investigate whether the system of equations derived above for z(s) and D(s) 
can lead to a distance-redshift relation close to the one observed. For wall universes we 
consider, 

For a chosen density contrast h(x) we can determine b(t, x) and a(t, x) and solve the photon 



geodesic Eq. (3.5) for a given angle 6q of the observed photon w.r.t. the y-plane, 



COS O = _^L±J_ (3.15) 

&(x ,t )n°(0) V ' 



We again set the initial value or the affine parameter to 0, hence xo = x(0) etc. 

2 

9t; 



We have investigated two choices for M(x). The first is simply M(x) = ^ (1 + e cos(fcx)) 



n 

which we have already discussed before. The results for this case are shown in Fig. |2[ 

The result is quite striking: The deviation from the Einstein-de Sitter distance-redshift 
relation is very small. On the level of a few percent in the most extreme case. Much smaller 
than the deviation for an open (Milne) Universe or even for ACDM. Hence voids and walls 
with the chosen parameters cannot simulate the observed distance redshift relation. We 
have also studied different values of the parameters (e, k), but all cases which are such that 
there is no singularity before to l ea d to small deviation from Einstein-de Sitter. Only for 
wavelengths of approximately Hubble scale, k ~ Hq, where we can choose e ~ 1CP 3 do the 
deviations become relatively large. But the density profile chosen here does not at all lead 
to a relation that resembles the observations. 

As a second profile we consider thin, highly concentrated over-dense walls with an 
exponential profile: 

hW =7ib£ exp ( z ^)- 1 ' (3 - 16) 
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Figure 2: We show the relative luminosity distance redshift relation ^ u ^ z ,\ = ^ n L>E . dS . {z) for 
different models with luminosity distance D(z). The blue dotted curve is for a Milne Universe, the 
red dashed curve is for ACDM universe with J1a = 0.7 and CIm = 0.3. The remaining two lines 
are our wall universe. The black solid line is in an under density while the purple dot-dashed line 
is in an over density. In the top panel, we consider light propagating in the a;-direction only. The 
bottom panel is the same but for light propagating in the y-direction. The parameters for the wall 
model are the physical ones, e = 9.5 x 10 -6 and A = 80 Mpc. 



VJ.U 

0.4 
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where A = Xj+i — Xj. In the limit a <C A the mean of h (x) vanishes and min x h (x) = — 1. 
Again, we choose e such that there is no singularity before to. The results for this profile 
are shown in Fig. ||. 

We have obtained the following result in these two examples (and other profiles which 
we do not present here explicitly): The modification of D(z) never goes beyond the case of 
the open universe. We do not obtain acceleration by a series of dense walls. Even though we 
present here only two simple profiles, we think the conclusion is valid beyond these cases: 
if a photon passes through many compensated under- and over-densities in the integrated 
distance D(z) the effect is minute as long as the time the photon spends inside a wall is 
much smaller than the time scale at which the gravitational potential of the wall evolves. 
A perturbative (first order) calculation gives a flavour of this effect. Indeed, at first order 
in the perturbed direction, the difference between D(z) in our models and DEds{ z ) OI a 
matter dominated universe can be written as 

d l ( Ze ) - Df dS = (i + z e ) fa - Ti e ) (| (h ( Vo ) + h ( Ve )j) - (i + z e ) j V ° jh (n) d v 

+ (1 + z e ) d V jf dn'^h" (rf) rf - p d V ^h" (rj) f}, (3.17) 

where the subscripts e and o respectively mean that the conformal time is evaluated at the 
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Figure 3: We show the relative difference between the distances in 'realistic' wall models and 
in EdS universe for photons propagating in x-direction. The top panel is obtained with e = 10 -9 , 
A = 40 Mpc and a = 1 Mpc, while the bottom one with e = 5x 10~ 8 , A = 15 Mpc and a = 1 
Mpc. In both cases the observer is at the center of the void. We have checked that the order of 
magnitude does not change for an observer in an over density. In the second case, we see that we 
obtain an effect of the same order of magnitude as the swiss cheese universe discussed in [E4| . 



source (emission) or at the observer and expresses the perturbation of the energy density 
in under and over densities (see Appendix C for a derivation of the linearized result). From 
this expression, valid in the linear regime only, and for a periodic perturbation, it becomes 
clear that the deviation of Dl(z) with respect to Df dS depends on the amplitude e of the 
perturbation and on the values of the conformal time at the source and at the observer. In 
the case of periodic perturbations, the contributions from photon path are mostly cancelled 
in the integral terms. Of course in the full non-linear calculation there is no simple relation 
between the matter over density h and the gravitational potential. In this case in principle 
the full non-linear Einstein equation have to be solved and Eqs. (|3.5| ) and ([O]) govern 
D L (z). 

Surprisingly, however, our non-linear simulations show that this result holds also to 
some extent in the non-linear regime. Note that, even though our value of e is small, the 
over densities in the walls are large at late times, such that they develop singularities soon 
after today and we are deeply in the non-linear regime. While we do not have a proof that 
our conclusion holds in all cases, we have tested this also with other periodic wall profiles. 

In Fig. ^ we show the deviations of the expansion rates with respect to the Hubble 
expansion in EdS universe. We note that the deviations in the unperturbed directions 
are small. However, in the perturbed direction these deviations can be large locally inside 
a wall, and they would be measurable by direct, local measurements of H{z). However, 
they compensate when averaged over a wall thickness and do not show up in integrated 
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quantities like D(z). 
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Figure 4: We show the relative differences between the expansion rates in the thin, highly 
concentrated over-dense wall model and the Hubble expansion in EdS universe. The top panels are 
obtained with e = 10~ 9 , A — 40 Mpc and a = 1 Mpc, while the bottom ones with e = 5 x 10~ 8 , 
A = 15 Mpc and a = 1 Mpc. In both cases the observer is at the center of the void. The left panels 
show the expansion rates in the perturbed direction, while the right ones in the y-direction. The 
results for the cosine profile not shown here are similar to the two top panels. 



3.3 Mimicking dark energy 



Yoo et al. [18] have shown that in an LTB model every given distance-redshift relation 
can be mimicked by a suitable choice of the density profile. The same is true for a wall 
universe. For a given function D(z) we can find a density profile which leads to exactly 
this distance-redshift relation for a photon coming in x-direction. First of all, for such 
a photon the shear vanishes for symmetry reasons and 1Z is given by (|1D. To find the 
density profile, which is equivalent to finding M(x) or M{z) = M(x{z)) we have to solve 
the following coupled system of six ordinary differential equations (in principle none of the 



other equations couples to ( 3.2C| ) since both, Fm and Fp do not depend on x explicitly), 
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which is very similar to the system solved in Ref. [18]: 



dM 

ds 
dp_ 

ds 
dx 

ds 
dt 

ds 
dz 

ds 

dC 

ds 



F M (t,z,M,P,0, 


(3.18) 


Ff,(t,z,M,p,0, 


(3.19) 


F M (t,z,M,p,Q 
P 


(3.20) 


l + z, 


(3.21) 


c 

dP ' 

dz 


(3.22) 


-4vr (1 + zf pD, 


(3.23) 



where we have defined 



c 



dz dD 
ds dz 



and f3 



dM 
ds 
dx 
da 



M' 



3*§ 



E. 



(3.24) 



In Appendix [A| we give the derivation of this system and the detailed expressions for Fm and 
F@. There, we also explain the method used to specify the initial conditions at the observer. 
All the constraints are fixed by requiring the system to have no critical points. Note also 
that z(s) need not to be monotonic. If dz/ds = at a value of s where £ = dD/ds ^ 0, the 
derivative dD/dz is not well defined. This is, however, not the case of a ACDM Universe 
which we want to mimic here. We are then left with one initial condition, which we choose 
by requiring 



(3.25) 



SO 



i.e. the value of the Hubble rate at the observer today does not depend on direction. 
In Fig. |5] we show M(x) as well as its derivative with respect to the x coordinate, j3 (x), 
for the solution mimicking the ACDM expression for D(z), for 0,k = 0, fl m = 0.3 and 
^pe(z) = 0.7 =constant. 



D(z) 
Xk{t) 
H(z) 



1 



l + z 

1 



XK 



dz' 



o H(z 
sm(r\^K) , and 

H (n m (i + zf + n K (i 



where 



(3.26) 



fi r (l + zf + n DE (z) 



1/2 



In Fig. ^, we show how the luminosity distance deviates when the observer looks at 
photons coming in with different angles #o- F° r $o = 90 degrees, we have photons traveling 
in x-direction, in this case the luminosity distance is fitted to the one of ACDM by solving 
the system of Eqs. (3.18- 3T23| ) with the functions M(x) and /3(x) shown in Fig. [|. It is 
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Figure 5: We show the function M(x), top panel, and its derivative /3(x), bottom panel. In 
principle, there is a entire family of functions M(x) parametrized by the initial value M(0) = M 
that we are free to choose (appendix A). Here, we present the solutions corresponding to Hq = 
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Figure 6: We show the relative differences between luminosity distances for photons traveling in 
the a>direction (perpendicular to the walls) and photons observed with an ang lc O (see Eq.(p^)). 
From the top to the bottom, we respectively have 9q = 75, 60, 45, 30, 15, 5 degrees. 
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interesting to remark that a given angle of 6q £ [0; 90] degrees at the observer corresponds 
to an angle at the emission 9 e > 8q. This is a consequence of the spacetime geometry 
induced by the walls: due to the clustering in direction x, corresponding to 9 = 90°, its 
expansion slows down in time. 



5r 1 i 1 i 1 i 1 i 1 i 




z 



Figure 7: We show the ratio of our density profile to the Einstein-de Sitter one as a function of 
the cosmological redshift. 

In Fig. 0, we present the density profile corrected by the isotropic expansion rate, 
(1 + z)~ s p(z)/po, po = p(z = 0), obtained for our model to mimic ACDM luminosity 
distance. Finally in Fig. |8[ we plot the expansion rates in the longitudinal and transverse 
directions, H a = a/a and Hb = b/b. It is interesting to estimate roughly the features of the 
under density needed to fit ACDM luminosity distance. For example, if one considers the 
highest redshift for which we have data from supernovae, at around z ~ 1.7. This roughly 
corresponds to a size ~ Hq 1 . (Of course we have another data point from the CMB. 
The angular size of the acoustic oscillations provides an excellent measure of the angular 
diameter distance to the last scattering surface, z ~ 1090. But this is not very relevant in 
our context as the Universe is to a good approximation matter dominated from z = 2 to 
z = 1090.) An under density of the size of the order of the Hubble distance is necessary 
to mimic ACDM with our walls. Moreover, we can also determine the ratio of the energy 
density normalized at the observer to the energy density in an Einstein-de Sitter model at 
z ~ 1.7 which is about 4. At high redshift, z > 10 the anisotropy is very small and the 
Universe is close to a Friedmann Universe with about 5 times the matter density obtained 
from local estimates. 

3.4 Redshift drift 

In the previous section we have fixed M(x) to reproduce the distance redshift relation of 
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Figure 8: We show the relative expansion rates in the transverse and longitudinal directions as 
functions of the cosmological redshift. We use the following notation : AH a ,b — H a ,b — Hf, where 
Hp is the expansion rate in an Einstein-de Sitter universe, and ff a ,6 are the expansion rates in the 
longitudinal and transverse directions, normalized to the their values at the observer. 



ACDM universe. Of course, having one free function to play with, namely M(x), we expect 
to be able to fit one function, in our case D{z). If we now proceed to another, independent 
observable, we shall most probably not fit it. We have done this by looking at the redshift 
drift, defined as the rate of change of the redshift of a co-moving source per unit of observer 
time. In a Friedmann Universe the redshift drift is simply 

dz _ z(t s + At,) - z(t s ) 

— = km — = H (l + z)- H(z), 3.27 

where H{z) = H{t s ) and Hq denote the Hubble parameter at the source position at time t s 
and at the observer at the moment to- We have computed the corresponding function (for 
light rays in x-direction) from our solution M(x). The general expression for the redshift 
drift of a wall Universe in x-direction is (see Appendix |B|), 



dz 
dt 




il + z) I \" T \(1 + Z >y 2 dz' 



= -(l + z )J* ^Gp -^)f(l + z') - 2 dz' . (3.28) 

Since we do not require Mq = as in LTB model, we can in principle have a positive 
redshift drift at low redshift; but we do not obtain this for our best fit profile M{x) with 
tB(x) = 0. The result is compared with ACDM in Fig. ||. 
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Clearly the redshift drift for the two cosmologies are very different. We do have a 
second function to play with, the bang time is (a;), so that we could probably fix this 
observable. This has been done for LTB models in |25[j . However, as it is shown there, 
models which have both, the same redshift distance relation and the same redshift drift as 
ACDM can be ruled out with a third observable, the kinematic Sunyaev-Zel'dovich effect 
which comes from the recession velocity of clusters. 



4. Conclusions 



We have studied the effect of matter perturbations on the luminosity distance in a model 
with planar symmetry described by the metric fl2.l|). Considering 'realistic' walls we find 
that the effect from density inhomogeneities is very small, it nearly averages out. It leads 
to fluctuations of the luminosity distance around the 'background' distance, but not to a 
significant global shift. Our results (Fig. |3|) show that these fluctuations are due to matter 
inhomogeneities at the source and the observer positions, without any relevant contribution 
from the integrated effects of light propagation, like in the linear approach ( |3.17 ). Hence we 
can not mimic acceleration with many dense walls which grow by gravitational instability. 
Since we consider pressure-less matter only, the amplitude of density fluctuations is limited 
by the presence of singularities. This is a limitation of the model. 

After having shown that 'realistic' wall models can not reproduce the observed distance- 



18 



redshift relation, we have determined the density profile which can mimic it. We have fixed 
the free function of our model, M(x), to mimic the luminosity (or angular) distance of the 
ACDM universe. We have shown that the observation of the redshift drift can distinguish 
between this model and ACDM. Abandoning the assumption of an uniform bang time we 
could arrange the second degree of freedom, to fit the redshift drift too. We have 

found that the redshift drift in our model can be positive at low redshift, contrarily to the 
LTB model p|. 



With our solution M(x) we can fit ACDM distance for photons coming in x-direction 
for positive x only. This preferred direction corresponds to the radial incoming direction 
for LTB model. The deviation from ACDM for photons coming from different angles is 
typically a few percent (see Fig. ^). 
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A. Derivation of the system of differential equations and initial conditions 
A.l The system 



Here we derive in more detail the system (3.18) to ( p.23 ) and give the initial conditions 
used for the solution. 

Since we choose the photon affine parameter such that n°|o = 1 we have 

l + z(s)=n°(s) = ^. 

as 

Furthermore, the null condition for a light ray in x-direction implies 

dt \ 2 ( b' ' \ 2 ( dx \ 2 / 26' \ ^ / \ 2 

d~s) = \EJ \~ds~) = \3t%M J ) {lb) ' (A,1) 

The geodesic equation gives 

dz d t cl , ,o b , ,o , . _ - 

S = ^ = -a (1 + Z » =-i7 (1 + 3) • (A - 2) 

Hence, when the expansion in x-direction changes into contraction, a = 0, also dz/ds passes 
through zero. However, this does not happen in our case which mimics ACDM. Noting 
that geodesies in x-direction have no shear, the Sachs focusing equation yields 



d 2 zdD /dz\ 2 d 2 D 



. ^1 ^- T = -^Gp(l + z) 2 D, (A.3) 
as z dz \ds J dz z 

where we have used 1Z = AttG (1 + z) 2 p. We can now rewrite these equations in terms of 



the system (3.18) to (3.23). To find the functions Fm and Fa we first derive the following 
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useful relations 



t = Jt r, (A.4) 



d§_ 

T ' = r I 2-^- --—), (A.5) 
i ^ 3M / ' 



^-^ (5 - 2TS/) -^f 5 ' (A - 6) 
* = (S- *rS> - 2r>S») - f^|^ (2,' + . (A.7) 

Here 5' always indicates the derivative of S with respect to its argument r while as for all 
other functions of (t,x) the prime denotes the partial derivative w.r.t. x and the dot the 
one w.r.t. t. The null condition for the light ray can be written as 



with 



dM A dB 

—A 1 + -j-B 1 = ±1, 

ds ds 



The geodesic equation takes the form 



where 



d -^A 2 + f B2 = 0, (A.10) 
ds ds 



A, = (5 - 2rS>) - (1 + (* - *rS' - 2r>S») , (A.ll) 

C 16 2 M 2 32 Mr 2 

^ = "Igr^ ^ 5 " C 1 + *) 9^4 ^ + T5 ) ' (AJ2) 



with C = ^ = From this we infer 

FM(t,z,M M ) = ± AiB2 B 2 MBi , (A.13) 

W.>;M,e,0=T AiB ^ A2B{ (A.14) 

Since r is a function of M, j3 and t, we now have expressed everything in terms of our 
variables (t, z, M, (3,() and the given function D(z). Explicitly, Fm and Fp are given by 

3t% /6M\ 2 ' 3 if ifi + t 1 + < 2S ' + tS "> 
Fm = ± ^, _ - Sy + rSg ,_ Tga • (A.15) 



± 



4 

1 1 



18$ M 



^ 6M y/ 3 ^ff 1 + (1 + *) (f ~ rS' - 2t*S") 

\T) ~ ss' + tSS"-tS* • (A,16) 
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A. 2 Initial conditions 



Let us now turn to the initial conditions at so = 0. Without loss of generality we can set 
x(0) = 0. Clearly also z(0) = 0. From definition ( |3.24j ) we have 



as 



(A.17) 



s=0 



Since this is an initial condition for the Sachs focusing equation, we have consistently with 
our affine parameter normalization |]23|, p7|, 



C(o) 



i. 



From (3.22) we note that our system of coupled differential equations 
z cr defined by 

Z — Zcr 

For our ACDM parameters z cr ~ 1.6. To obtain a regular solution 
impose C, (z cr ) = 0. We remark that Eqs. ( A.1S| ) and ( |A.2j ) imply 



a - TT 

a 



(A.18) 

has a critical point 
(A.19) 

we must therefore 
(A.20) 



where we have used 



dD 

dz 



Hence the rate expansion in x-direction coincides with the measured 
In order to solve the system of five differential equations (Eq. ( |3.20| ) 



equation, since the solution x (s) can also be inferred from Eq. ( 3.21 ) via 
five initial conditions are needed. However, we only have two of them 



*(o) = c(o) 



1. 



Hubble expansion, 
is an independent 
the null condition), 

(A.21) 



We have two other constraints which we must satisfy at the critical point where ( A.1S| ) 
holds. Denoting the affine parameter at the critical point by s cr , we have 



C M = o. 



(A.22) 



These lead to two other initial conditions which can be determined using the shooting 
method. One remaining constraint is needed and we fix it by requiring 



Hn 



(A.23) 



This last condition fixes M(0) and makes sure that the Hubble rate measured today is the 
same in any direction. We then numerically integrate the system from the critical point to 
the observer by varying the three remaining conditions at the critical point until the initial 
conditions ( A.21 ) and ( A.23[) are satisfied. This matching is obtained by using the three 
dimensional Newton-Raphson method. Once the desired precision has been reached, the 
two remaining initial conditions /3(0) and i(0) can simply be read from the numerical data. 



21 



B. Derivation of the system of differential equations for the redshift drift 



The redshift drift for a LTB model has been derived in ||26|| . This approach can also be 
applied to our model. The null condition for the light ray (in ^-direction) and the geodesic 
equation lead to 

dz b . dt b' i — , 

We consider two infinitesimally close geodesies at fixed comoving position x, parametrized 
by 

{z c , t c } and {z c + 5z, t c + 5z] . 

Since the geodesic {z c ,t c } satisfies ( |B. 1|) , it follows 

d5z b' , . . b' 
-d^=E {l + z)5t+ E^ 

dst y . 



d5z b' 5z . 

dz b' I + z 

dSt 6t (B.3) 



Then, inserting (|B.l| ) we obtain 



Integrating ( |B.3|) we find 



dz 1 + z 



l + z 

This solution together with (B.2) leads to 



d ( 5z\ 1 (b' + 5z 



dz \5t J 1 + z \ b' ^ 



This equation is solved by ( p.28| ). Deriving the Einstein equation ( |2.3| ) twice (once w.r.t. 
x and once w.r.t. t), we obtain 

r, 2Mb' M' 
With (|J) and Q this results in the second line of (gjH). 



C. The linearized approach 

We determine the luminosity distance within linear perturbation theory for small deviations 
from a Friedmann-Lemaitre background. Let us define 

a(t,x) = a(t)(l + ef(t,x)), (C.l) 
b(t,x) = a(t)(l + eg(t,x)), (C.2) 
p(t,x) = p(t)(l + e6(t,x)), (C.3) 
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where the unperturbed quantities a (t), p(t) satisfy the Einstein equations for a flat matter 
dominated Friedmann universe (EdS). The perturbed quantities are determined by the 
Einstein equations at first order in e, 



;!/|/:; =87rGpS, (C.4) 

g' = o, (C.5) 

t 1 / 3 (2g + tg) = 0, (C.6) 

tfg" - tV* (2/ + 25 + 1 ( f + g) ) = 0. (C.7) 



Neglecting the decaying mode and imposing that at the beginning the scale factors in all 
three directions agree, we obtain [jDJ, 

9 = |, (C8) 

/ = ^o4 /3 * 2/3 + y, (C9) 

where 5o (x) = 5 (t, x) + / (t, x) + 2g (t, x) is independent of time. This is a consequence of 
energy conservation and can also be derived by combining ( |C.4| ) to ( |C . 7| ) . 

We are interested in finding the relation between 5o and M, E in the perturbative 
regime. Following [16| we expand the solution (2.7, |2.8|) around rj = in terms of tB ^ x ' <g 1 



and <C 1. Comparing the expanded solution with the linear one we find 



eS 

(1 + ed ) , V = 

With the ansatz (|2~T7| , gl9|) we can identify (5o (#) with /i (x) in the perturbative regime. 



M = -^(l + e5 ), E = ^. (CIO) 



The angular distance is determined by Sachs focusing equation (|3.6| ) . We note that the 
shear term does not contribute to first order. Since light propagation is not affected by a 
conformal transformation, it is convenient to work with the conformally related geometry 

ds 2 = -di] 2 + (1 + 2ef) dx 2 + (1 + leg) {dy\ + dyj) . (C.ll) 

From this, we compute the Christoffel symbols (here we denote the derivative w.r.t. the 
conformal time rj by a dot) 



r{ = ef, 



r?! = ef, 


1 22 


"pO ~ 

— 1 33 — 


eg, 


rli = ef, 


r 1 - 
1 22 - 


"pi rvj 

" 1 33 — " 


-eg' 


20 — 1 30 — e i/' 


r 2 
1 21 


"p3 ~ 

— 1 31 ~ 


eg', 



and the Ricci tensor 



#00 = -e(/ + 2$) , 
R 10 = -leg', 



R xl ^e[f-2g"), 
R22 = R 3 3 = e{g- g") 
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At O-order we are free to parametrize the affine parameter s such that n° = 1 and n l = 5 l1 
(we are interested in the distance in ^-direction) . With this we obtain the coefficient 1Z 

K = -e(g + g" + 2g>). 

Consistently with the parametrization of the affine parameter s such that n° (so) = 1, the 
initial conditions are D (s G ) = and D' (s ) = — 1. After an integration by parts we find 
the solution to Sachs focusing equation (3.6), 

D (s) = (s - s) (1 + eg (s ) + eg (s)) + 2 f ds'eg (s') . (C.12) 

With the above initial conditions for the Sachs focusing equation, we consider a thin light 
bundle with the vertex at the observer position. Hence the solution (C.12) is the angular 
diameter distance, see [23]. To determine the luminosity distance we have to compute also 
the redshift, using the geodesic equation for n , 

1 + z= 9^u»\ =1 _n dse ^ (ai3) 

where | e denotes the emission point, the source, and we denote the affine parameter at 
the source by s e . With the same geodesic equation we derive the relation between the 
conformal time rj and the affine parameter s, n° = drj/ds, 



rio-Ve = So- s e + ds ds'ef(s'). (C.14) 

J S J S Q 

In terms of conformal time the luminosity distance then becomes 

Dl (Ve) = {Vo ~ Ve) [ 1 + eg D + eg e - 2 / d V ef ) +2 / dr] eg - di] dr/'ef. 

V J Vo / J Vo J Vo J Vo 

(C.15) 

All of this is valid in the conformal geometry, where the expansion of the Universe is divided 
out. Taking into account the expansion of the universe, changes the relation between the 
affine parameter and conformal time. The luminosity distance scales as || 

Dl = —, , D L = —— = (1 + z e ) D L . 
a (rj e ) a (rj e ) 

Since conformal time is not an observable quantity, we rewrite the distance in term of the 
observed redshift. We define the observed redshift as z e = z e + 5z e and we compute the 
correction term. The same calculation as presented in Ref. || leads to 

^D-L^j 5z e = (( Vo - rie) + H' 1 ) 5z e , (C.16) 

where 

5z e = -(l + z e ) d7]ef. (C.17) 

J Vo 
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Subtracting QC.16D we obtain the distance-redshift relation 

D L (z e ) = (1 + z e ) (rj - r) e ) (l + e 9o + eg e - J* dr, ef\ (C.18) 

( r Ve r Ve p A i + 2! f Ve 

+ (1 + ^) 2/ drjeg- dr, drfef + dr, ef. 

V J rio Jrio Jri J "e Jri 



With 



/•»7e rv . /•'yo . r»?o rv 

- dr, drfef = (r, e - r, Q ) I dr, ef + dr, dr,' 'ef ', 

Jr] Jr\ Jr) e J r) e J r) e 

we can rewrite the above expression in the form as 

D L (z e ) = (1 + z e ) (r, a - r, e ) (1 + eg Q + eg e ) (C.19) 
/ rvo rvo rv A i + z rV» 

+ (1 + z e ) -2 / dr, eg + / dr? / drfef - / dry ef . 



Using the solutions ( |C.8| , |C.9| ) we express the distance in terms of 5o (r,). Conformal time 
is defined as 

dr, = => »/(*) = 3t 1/3 to /3 , setting 77 (0) = 0. 

This leads to 

/(r?,x(r / )) = l^(x(r ? ))r ? 2 + ^M, 
f(r,,x(r,)) = (x(r,))r,, 
and consequently to the following distance-redshift relation 



D L (z e ) = (1 + z e ) (r,o ~ Ve) (l + | (<*o (%(Vo)) + $o (x(Ve)))) - (1 + z e ) J j$o (x(rj)) dr, 



+ (1 + z e ) I" dr, f" drfUl (x(rf)) rf - ±±*L P drj (x(r/) ) v . (C.20) 

References 

[1] A. Riess ct al. Astron. J. 116, 1009 (1998) [arXiv:astro-ph/9805201vl]; 

S. Perlmutter et al, Astrophys. J. 517, 565 (1999) [arXiv:astro-ph/9812133vl]. 

[2] R. Amanullah et al., Astrophys. J. 716, 712 (2010) [arXiv:1004.1711vl]; 
J. Guy et al. (2010), arXiv:1010.4743vl; 

T. Holsclaw, et al, Phys. Rev. Lett. 105, 241302 (2010) [arXiv:1011.3079vl]; 
N. Suzuki et al. (2011), arXiv:1105.3470vl. 

[3] J. Dunkley et al., Astrophys. J. 739, 52 (2011) [arXiv: 1009.0866]. 



25 



[4] R. Durrer, Phil. Trans. R. Soc. A 369, 5102 (2011), doi: 10.1098/rsta.2011.0285 
[arXiv:1103.5331]. 

[5] C. Blake et al., Mon. Not. Roy. Ast. Soc. 418, 1725 (2011) [arXiv: 1108.2637]. 

[6] C. Bonvin, R. Durrer and A. Gasparini, Phy. Rev. D 73, 023523 (2006) [arXiv:astro- 
ph/0511183v4]. 

[7] C. Vale and M. White, Astrophys. J. 592, 699 (2003) [arXiv:astro-ph/0303555vl]; 
K. Kainulaincn and V. Marra, Phys. Rev. D83, 023009 (2011) [arXiv: 101 1.0732]. 

[8] S. Rasanen, Class. Quant. Gravity 28, 164008 (2011) [arXiv: 1102.0408]. 

[9] G. Lemaitrc, 1933 Ann. Soc. Sci. Bruxelles A 53, 51 (1933); (English translation, with histor- 
ical comments: Gen. Rel. Grav. 29 637, 1997). 

[10] R.C. Tolman, Proc. Nat. Acad. Sci. USA 20, 169 (1934); (reprinted, with historical comments: 
Gen. Rel. Grav. 29 931, 1997). 

[11] For recent overviews see: 

K. Bolejko, M-N. Celerier and A. Krasinski, Class. Quantum Grav. 28, 164002 (2011) 
[arXiv:1102.1449v2]; 

V. Marra and A. Notari, Class. Quantum Grav. 28, 164004 (2011) [arXiv:1102.1015]. 

[12] P. Szekeres, Com. Math. Phys. 41, 55 (1975); Phys. Rev. D12, 2941 (1975). 

[13] A. Krasinski and K. Bolejko, Phys. Rev. D 83, 083503 (2011), [arXiv:1007.2083v2]. 

[14] N. Meures and M. Bruni, MNRAS 419, 1937 (2012) [arXiv:1107.4433v2] (DOI: 10.1111/j.l365- 
2966.2011.19850.x). 

[15] H. Collins (2010), arXiv:1011.2046vl. 

[16] A. V. Zakharov, Soviet Physics Journal 30, 1015 (1987). 

[17] J. Plebanski and A. Krasinski, An Introduction to General Relativity and Cosmology, (Cam- 
bridge, 2006). 

[18] C.-M. Yoo, T. Kai, K. Nakao, Prog. Theor. Phys. 120, 937 (2008) [arXiv:0807.0932]. 

[19] E. Komatsu et al., ApJS, 192, 18, (2011) [arXiv:1001.4538v3]. 

[20] H. El-Ad and T. Piran, ApJ, 491, 421(1997) [arXiv:astro-ph/9702135vl]. 

[21] P. Schneider, J. Ehlers and E.E. Falco, Gravitational Lenses, (Springer Verlag, 1992). 

[22] N. Straumann, General Relativity with Applications to Astrophysics, (Springer Verlag, 2004). 

[23] S. Seitz, P. Schneider and J. Ehlers, Class. Quantum Grav. 11, 2345 (1994) [arXiv:astro- 
ph/9403056vl]. 

[24] V. Marra, E. W. Kolb, S. Matarrese and A. Riotto, Phys. Rev. D 76, 123004 (2007) 
[arXiv:0708.3622v3]. 

[25] Philip Bull, Timothy Clifton, Pedro G. Ferreira, arXiv:1108.2222 (2011). 

[26] C.-M. Yoo, T. Kai, K. Nakao, Phys. Rev. D 83, 043527 (2011) [arXiv: 1010.0091]. 

[27] V. Pcrlick, Living Reviews in Relativity, 7, 9 (2004). 



26 



